Size Effects in the Magnetoresistance of Graphite: Absence of Magnetoresistance in 

Micrometer size Samples 
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We present a study of the magnetoresistance of highly oriented pyrolytic graphite (HOPG) as 
a function of the sample size. Our results show unequivocally that the magnetoresistance reduces 
with the sample size even for samples of hundreds of micrometers size. This sample size effect is 
due the large mean free path and Fermi wavelength of carriers in graphite and may explain the 
observed practically absence of magnetoresistance in micrometer confined small graphene samples 
where quantum effects should be at hand. These were not taken into account in the literature yet 
and ask for a revision of experimental and theoretical work on graphite. 
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Graphitic systems are nowadays a field of intensive 
activity There have been observations of quan- 

tum Hall effect in HOPG 0, Q as well as very large 
anisotropy in the electrical conductance (ratio between 
current parallel to perpendicular to the graphene planes) 
larger than 10 4 at room temperature Graphite looks 
as a good conductor in plane and an insulator between 
planes leading to a weak screening to external electric 
fields 0]. This means that an external electric field 
penetrates by tens of nanometers in graphite, in con- 
trast to a normal metal where the held is screened in 
the hrst atomic layers. In fact the dielectric constant 
of graphite at optical frequencies is positive, insulator 
like (e planc = 5.6 + i7.0,e c = 2.25) @, @. Electric 
held microscopy detects that regions of graphite are 
more insulating than others upon the interconnections 
of graphite planes produced by defects and their over- 
all density that influences the density of states and the 
Fermi level. The material exhibits a huge magnetic field 
driven metal-insulator transition [l], @|- Its band struc- 
ture and intcrband transitions with electron and hole 
carriers manifest in a huge ordinary magnetoresistance 
(OMR)@, H|. All these effects happen m macroscopic 
size samples of the order of millimeters. 

A large research activity has been recently started on a 
few graphene layers (FLG) [l(| with typical size of a few 
microns. Strikingly, the OMR in FLG samples is prac- 
tically suppressed even at T < 4 K, in contrast to bulk 
HOPG or Kish graphite where the OMR is > 1000% at 
B > 0.5 T at low temperatures Earlier work with 
graphite samples in the micrometer range showed sim- 
ilar behavior 
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[12j. However, the question on what 
happens when the graphite sample size is reduced has 
not been correctly addressed and the experimental data 
may need a new interpretation. The size of the sample 
is very important for defining the properties of the sys- 
tem because the de-Broglie wavelength for massless Dirac 



fermions Ad ~ hvp/Ep (vp — 10 6 m/s is the Fermi ve- 
locity and a typical Fermi energy ksEp < 100 K) or 
for massive carriers with effective mass m* < O.OLn, 
A m = h/\/2m*Ep, as well as the Fermi wavelength 
\p ~ 27m" 1 / 2 are of the order of microns or larger due to 
the low density of Dirac and massive fermions. Moreover, 
one can show that in such micron size systems new quan- 
tum mechanical oscillations appear as hne structure su- 
perimposed to the usual Schubnikov de-Haas (SdH) mag- 
netoresistive oscillations O, an experimentally observed 



but not yet recognized fact [1J, [15J, [16|, [ljj ■ The aim of 



this paper is to analyze experimentally the behavior of 
the OMR in HOPG samples upon their macroscopic size 
and the effects of a constrained region. We found that 
the smaller the size of a sample or constraint the smaller 
is the OMR. Surprisingly, the effect is already noticeable 
at large sample sizes, hundred of microns. 

In order to carry out a systematic study we have per- 
formed experiments in different HOPG samples from Ad- 
vanced Ceramics or Structure Probe Inc. with mosaici- 
ties 0.4° ± 0.1°, 0.8° ± 0.15° and 3.5° ± 0.4°. The smaller 
the mosaicity the more ordered is the sample. The re- 
sults presented in this letter are from HOPG samples 
with mosaicities of 0.4° and 0.8°. Results for larger mo- 
saicity will be presented elsewhere. The advantage of us- 
ing HOPG of good quality is that due to the perfection of 
the graphene layers and low coupling between them, a low 
two-dimensional carrier density n ~ 10 10 . . . 10 11 cm" 2 is 
obtained 0, EH • This value is much smaller than in typ- 
ical FLG's probably due to lattice defects generated by 
the used method to produce FLG samples and/or surface 
doping. Other advantage of using HOPG samples is that 
the preparation procedure for macroscopic sample size 
(from mm to (j,m) is rather simple enhancing the repro- 
ducibility and also the possibilities of checking different 
geometries. 

The graphite samples of thickness between ~ 10 /im 
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and ~ 40 /im were prepared by peeling from HOPG 
using double sided tape. The samples were then glued 
with varnish onto silicon substrates. Electrical contacts 
were made by using silver loaded cryogenic epoxy in 
the usual four-point arrangement for resistance measure- 
ments. The current applied was between 1 /iA to 100 mA. 
The sample width was reduced by cutting the sample 
with a diamond tip (wide cuts). Also we changed the 
voltage electrode distance without changing the sample 
size. For small constrictions we used a FEI Novalab 200 
dual beam scanning electron microscope (SEM) with a 
crossing focus ion beam (FIB). This set up allows us de- 
fined production of rectangular constraints of different 
lengths and widths between the voltage electrodes and 
studied their influence on the OMR. We have produced 
constrictions of micrometer size and not in the nanometer 
range because we are interested to know at what sample 
size or constriction width one starts observing size effects. 
Part of the experimental studies were done at room tem- 
perature and 77 K and a few samples were cooled down 
to 4 K. The magnetic field was applied always perpen- 
dicular to the graphene planes of the HOPG samples. 

Electron backscattering diffraction (EBSD) was per- 
formed with a commercially available device (Amctek- 
TSL). In this set up the sample under investigation was 
illuminated by the SEM beam and the diffracted elec- 
trons were detected by a fluorescence screen and a digi- 
tal camera. The TSL software was used to calculate the 
orientation of the HOPG surface as function of the elec- 
tron beam position. Figures 1(a) and (b) show the same 
600 x 200 /im 2 area of a HOPG sample. Figure 1(b) 
shows the distribution of the c-axis of the hexagonal unit 
cell at the graphite surface depicted by the spreading of 
the (red) color. The grain distribution can be also seen 
in Fig. 1(a) where the in-plane orientation is recognized 
by the (blue-green) color distribution. In Fig. 1 we see 
that the typical crystal size in HOPG (0.4°) is of the or- 
der of micrometers. Taking into account this fact and 
that upon defect concentration in the crystals and their 
influence on the effective carrier density Xd,™. ^ 1 /^rn, 
we expect that size effects might be seen at relatively 
large sample or constriction size. This should be in prin- 
ciple possible also because the electron mean free path 
L e ~ m*VF^/e ^ 1 M m m HOPG, taking into account 
a mobility u~ 10 6 cm 2 /Vs measured in samples of low 
mosaicity |l8|, a value much larger than those found in 
typical FLG. All these numbers already suggest a sit- 
uation very different from that found in metals, where 
size effects in the magnetoresistance may start to be seen 
when the sample size reduces to ~ 20 nm or less. 

Figure [2] shows the dependence of the OMR defined 
as (R(0.5T)-R(0))/R(0) with the sample width mea- 
sured at 300 K (•,□) and 77 K (*) (right axis). Note 
that the sample width changes between 1 cm and 80 /im. 
The OMR clearly decreases from ~ 40% for a macro- 
scopic sample to ~ 20% for a still relatively large con- 




FIG. 1: Inverse pole figures maps from EBSD measurements 
of the HOPG surface showing the in-plane (a) and the c-axis 
orientation (b) of the crystallites on the sample surface. The 
scan size is ~ 600 x 200 /im 2 ; the size of the domain patterns 
matches the results of EFM measurements 0]. 



stricted sample of 80 /im. The data (•) also indicate 
that the OMR tends to saturate for a sample width above 
~ 300 /im. In Fig. [5] there are points for another HOPG 
sample (□) prepared similarly, which shows the same 
trend. The decrease of OMR is also clearly observed 
at 77 K (*). This result is quite spectacular because the 
reduction is a factor of 2 at 300 K and nearly a factor 10 
at 77 K, even if the sizes are ~ 50 times the (estimated) 
mean free path or the Fermi wavelength. One may spec- 
ulate that due to the method used to reduce the width 
of the sample we are introducing defects and this is the 
reason for the observed decrease of the OMR. To check 
that this is not the case we have done two more exper- 
iments that are described below. The reason is that by 
reducing temperature the mean free path increases and 
quantum, non-classical effects start to affect larger areas 
of the sample. Note that graphite has Dirac and massive 
electrons and the former have a tremendously large mean 
free path (several tens of microns) and these are likely the 
ones that reduce the OMR at large sample size. 

In the same Fig. [2] we show the change of OMR at 
300 K for a HOPG sample changing only the distance 
between voltage electrodes from ~ 1 cm to ~ 100 /im 
(~k). The same tendency is observed, i.e. the OMR de- 
creases nearly logarithmical with the sample width or 
electrode distance. This phenomenon is also observed 
if we produce a rectangular constraint of small length 
~ 2 /im, much smaller than the sample width ~ 500 /im, 
and located at the middle between the voltage electrodes 
(distance ~ 500 /im) on a HOPG sample (0.4°). The 
constraint length is so small that in principle one does 
not expect any change in the OMR. In fact, at 270 K the 
value obtained for the OMR for this sample is ~ 0.30, in 
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FIG. 2: Magnetoresistance as a function of sample's width 
size in micrometers. The data were taken from three HOPG 
samples of mosaicity 0.4° at 300 K (•,D,-^ r ), left axis, and 
a sample of 0.8° at 77 K (*), right axis. The width of the 
samples described by the points (•,□,*) were changed using a 
diamond tip. The straight line is only a guide to the eye. The 
width values of the sample described by the symbol (if) mean 
the distance between the voltage electrodes without changing 
the sample dimensions. Inset: Magnetoresistance R(B)/R(0) 
vs. applied field for a HOPG sample with a constraint width 
of 16 /im at 270 K and 4 K. The dashed (red) line follows the 
function 70IBI 1 ' 4 + 1. 



agreement with the data shown in Fig. [21 and remains 
independent of the width of the narrow-constraint. In 
other words at this temperature we measure the OMR 
coming from the sample bulk. The inset in Fig. [5] shows 
the parabolic behavior of the OMR at 270 K for the sam- 
ple with a constraint-width of 16 /im; this relationship 
remains independent of the constraint-width in the mea- 
sured range. 

However, at lower temperatures, due to the increase of 
the mean free path the influence of the constraint might 
be not negligible anymore. In agreement with this expec- 
tation we observed a clear dependence of the OMR with 
constraint width at 4 K, see Fig. [3) The inset in this fig- 
ure shows the magnetoresistance defined as the resistance 
ratio i?(0.5)/i?(0) vs. the width of the constraint where 
a logarithmic dependence is clearly observable. Two im- 
portant details of these results are worth to mention. 
The observed influence of the constraint remains even for 
macroscopic large widths, see inset in Fig. [3] Note that 
the measured magnetoresistance is given by two contri- 
butions, one coming from the sample bulk, which remains 
as a constant background, and the other due to the influ- 
ence of the constraint. This explains the still relatively 
large OMR observed for constraint width of the order of 
~ 2 fim. 

All the results presented above point out to phenom- 
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FIG. 3: Magnetoresistance of a HOPG sample (0.4°) as a 
function of applied magnetic field at different widths of a 
constraint localized in the middle of the voltage electrodes, 
at 4K. The constraint widths were (1.8, 16, 35, 500) fim corre- 
sponding to the symbols: (A, -k, O, •)• The straight line has 
a field dependence B 12 and the dashed line B 1 ' 7 . Inset: Mag- 
netoresistance defined as R(0.5)/ R(0) vs. constriction width. 
The line is a logarithmic fit to the data. 



ena that in a metal do not happen at least at the scale 
of the graphite samples. Our results show that reduc- 
ing the sample size the OMR reduces. If we extrapolate 
the overall sample size to a few micrometers or below 
we expect a negligible OMR, as has been experimen- 
tally observed for graphite samples of size of the order 
of ~ 50 /im to less than 100 nm [TllE^]. Classical theory 
may not be applicable when the coherence length, mean 
free path, and the Fermi wavelength of the carriers are of 
the order of microns, the sample size. In these conditions 
Boltzmanns transport theory may not be a good approx- 
imation and all kind of finite size effects should appear. 
Strictly speaking we should ask: What is the meaning of 
the classical cyclotron radius r c ~ m*vp/eB ^0.1 ^m at 
B ~ 1 T when the wavelength of the electrons > 1 /im? 
The Onsager quantization that gives rise to the de-Haas- 
van Halphen and SdH oscillations needs to be also re- 
considered because now the wavevector might not be a 
continuum, etc. 

The reported OMR data in FLG indicate practically 
the absence of OMR [HEEIIZl 

in comparison with that 
observed in bulk HOPG, which is of the order of 5000% 
at ~ 0.5 T at low temperatures. The sample size in 
the experiments with FLG's is what strongly reduces the 
OMR and not their thickness. In fact, we have done 
experiments with HOPG samples of different thicknesses, 
from 1 mm to ~ 10 fim. Our results show clearly that 
reducing the thickness the OMR increases. 

There is still a behavior that we would like to remark 
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about the OMR. At 270 K the OMR increases quadrati- 
cally with field (see inset in Fig. [2]) while at temperatures 
below ~ 100 K it behaves quasi linearly (see, for example, 
Fig. [3]) . This is an experimental fact that is well know for 
graphite and other semimetals. Abrikosov TfJ 2(j sug- 
gested that this phenomena takes place in the Landau 
level quantization regime of the Dirac fermion system, 
i.e. above a certain minimum field. However, the fact 
that this quasi linear OMR is also observed at very low 
fields (< 0.01 T, see for example the data (■) in Fig. [3j 
casts doubts about this interpretation. Our explanation 
is again in term of electron coherence and the role played 
by the platelets of ~ 10 fim size approximately conform- 
ing HOPG (sec Fig.l). When the mean free path of the 
electrons is larger than the platelets then electrons re- 
flect and are being transmitted coherently at the borders 
of the platelets and quantum currents circulate creating 
Hall-like potentials. This circulation of quantum cur- 
rents at the platelets edges (see Fig. 1) produce a quasi 
linear behaviour with field of the OMR. A similar pic- 
ture has been discussed to interpret experiments done 
with a Corbino disk geometry in HOPG [2l| and disor- 
dered semiconductors (22|. This is the same thing that 
happens in optics with the reflection of waves between 
different media with the addition that some Dirac elec- 
tron existing in graphite are dispcrsionlcss as photons. 
At room temperature the mean free path is smaller than 
the platelets size and the quantum currents at the border 
of the platelets do no exist. In this case we recover the 
classical parabolic ohmic behavior of the OMR. This ap- 
pears to be the explanation for the two regimes of OMR 
at low and high T. The OMR 4 K-data for different con- 
straint widths shown in Fig. [3] indicate that the low-field 
slope increases (from 1.2 to 1.7) decreasing the size of the 
constraint. This behavior is due to electron localization 
at small fields that shows a OMR dependence of the type 
<~ B~ 5 . The change of slope in Fig. [3] is a clear mani- 
festation of this effect. Localization effects are reflected 
also in a noise behavior of the OMR at low fields. 

In conclusion, our studies reveal a rather spectacular 
size effect in the transport properties of graphite, spe- 
cially in its magnetoresistance. The influence is observ- 
able up to macroscopic sample size of the order of several 
hundreds of micrometers. We suggest that this is due to 
the extremely large coherent and Fermi wavelength of the 
carriers (Dirac as well as massive electrons) in graphite. 
It seems clear that when the sample size is of the order 
of the Fermi wavelength and smaller than the carriers 
mean free path the classical Boltzmann transport theory 
becomes gradually not valid. In particular for these sizes 
the wavevector of the carriers might conform a discrete 
set of quantum numbers and therefore the quantization 
of orbits requires to have this into account. Our observa- 



tions provide a possible answer to the practically absence 
of OMR in FLG's reported in the literature. 

This work was done with the support of the DFG under 
ES 86/11 and of the Spanish CAICyT. One of us (J.B-Q.) 
is supported by the EU. 



* Electronic address: 
' Electronic address: 



[1 

[2 

[3 

[4 

[5 
[6 

[< 

[« 

[9 

[10 
[11 

[12 

[13 

[14 

[15 



nicolas . garciaQfsp . csic.es 



esquin@physik . uni-leipzig. de 



Y. Kopelevich, P. Esquinazi, J. H. S. Torres, R. R. 
da Silva, and H. Kempa, Advances in Solid State Physics 
43, 207 (2003). 

Y. Kopelevich, J. H. S. Torres, R. R. da Silva, F. Mrowka, 
H. Kempa, and P. Esquinazi, Phys. Rev. Lett. 90, 156402 
(2003). 

R. Ocana, P. Esquinazi, H. Kempa, J. Torres, and 

Y. Kopelevich, Phys. Rev. B 68, 165408 (2003). 

J. Gonzalez, F. Guinea, and M. A. H. Vozmediano, Phys. 

Rev. B 63, 134421 (2001). 

H. Hoinkes, Rev. Mod. Phys. 52, 933 (1980). 

N. Garcia, W. E. Carlos, M. W. Cole, and V. Celli, Phys. 

Rev. B 21, 1636 (1980). 

Y. Lu, M. Munoz, C. S. Steplecaru, C. Hao, M. Bai, 
N. Garcia, K. Schindler, and P. Esquinazi, Phys. Rev. 
Lett. 97, 076805 (2006). 

H. Kempa, Y. Kopelevich, F. Mrowka, A. Setzer, J. H. S. 
Torres, R. Hohne, and P. Esquinazi, Solid State Com- 
mun. 115, 539 (2000). 

W. Zhi-Ming, X. Ding-Yu, Z. Shi- Yuan, X. Qing-Yu, 
M. VanBael, and D. You-Wei, Chin. Phys. Lett. 24, 199 
(2007). 

M. I. Katsnelson, Materialstoday 10, 20 (2007). 

E. Dujardin, T. Thio, H. Lezec, and T. W. Ebbesen, 

Appl. Phys. Lett. 79, 2474 (2001). 

X. Zhang, Q. Z. Xue, and D. D. Zhu, Phys. Lett. A 320, 
471 (2004). 

N. Garcia, |http:// arxiv.org/find/cond-mat/0706.0135 



to be published. 

I. A. Luk'yanchuk and Y. Kopelevich, Phys. Rev. Lett. 
93, 166402 (2004). 

K. S. Novoselov, A. K. Geim, S. V. Morozov, S. V. 
Dubonos, Y. Zhang, and D. Jiang, Science 306, 666 

(2004) . 

[16] K. S. Novoselov, A. K. Geim, S. V. Morozov, S. V. 
Dubonos, Y. Zhang, and D. Jiang, Nature 438, 197 

(2005) . 

[17] Y. Zhang, Y.-W. Tan, H. Stormer, and P. Kim, Nature 

438, 201 (2005). 
[18] M. S. Sercheli, Y. Kopelevich, R. R. da Silva, J. Torres, 

and C. Rettori, Solid State Commun. 121, 579 (2002). 
[19] A. A. Abrikosov, Phys. Rev. B 60, 4231 (1999). 
[20] A. A. Abrikosov, Europhys. Lett. 49, 789 (2000). 
[21] H. Kempa, P. Esquinazi, and Y. Kopelevich, Solid State 

Communication 138, 118 (2006). 
[22] M. M. Parish and P. B. Littlewood, Nature 426, 162 

(2003). 



